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Abstract 

Gauge theory-string theory duality describes strongly coupled Af = 4 supersymmetric 
SU (n c ) Yang-Mills theory at finite temperature in terms of near extremal black 3-brane 
geometry in type IIB string theory. We use this correspondence to compute the leading 
correction in inverse 't Hooft coupling to the shear diffusion constant, bulk viscosity 
and the speed of sound in the large-n c M = 4 supersymmetric Yang-Mills theory 
plasma. The transport coefficients are extracted from the dispersion relation for the 
shear and the sound wave lowest quasinormal modes in the leading order a'-corrected 
black D3 brane geometry. We find the shear viscosity extracted from the shear diffusion 
constant to agree with result of |hep-th/0406264 ; also, the leading correction to bulk 
viscosity and the speed of sound vanishes. Our computation provides a highly nontrivial 
consistency check on the hydrodynamic description of the a'-corrected nonextremal 
black branes in string theory. 
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1 Introduction 



The correspondence between gauge theories and string theory of Maldacena [1,2] has 
become a valuable tool in analyzing near-equilibrium dynamics of strongly coupled 
gauge theory plasma [3-19]. The best studied example of strongly coupled thermal 
gauge theory plasma is that of the Af = 4 SU (n c ) supersymmetric Yang-Mills theory 
(SYM). In the large- n c limit, and at large 't Hooft coupling g\ M n c 3> 1, the holographic 
dual description of the M = 4 plasma is in terms of near-extremal black 3-brane 
geometry in type IIB supergravity [20]. In this case one finds [5,6,8] the speed of 
sound c s , the shear viscosity rj, and the bulk viscosity £ correspondingly 

Cs = i=, V = \nlT\ C = 0. (1.1) 

In a hydrodynamic approximation to near-equilibrium dynamics of hot gauge theory 
plasma there are several distinct ways to extract the transport coefficients (jl.ip . First 
[5], the shear viscosity can be computed from the two-point correlation function of the 
stress-energy 1 tensor at zero spatial momentum via the Kubo formula 



r, = lim -L / dtdx e iut ([T xy (x), T xy (0)]) . (1.2) 

w— >0 luj J 



Second [6], the diffusive channel two-point retarded correlation function of the stress 
energy tensor, for example, 

G tx , tx {u,q) = -i [ dtdxe^- i(lz e{t){[T tx {x), T ta (0)]> oc . * 2 (1.3) 

has a pole at 

uj = -iVq 2 , (1.4) 
where the shear diffusion constant T> is 

* = (1.5) 

with s being the entropy density of the gauge theory plasma. From the thermal field 
theory perspective it is clear that computation of the shear viscosity via Kubo rela- 
tion (jl.2j) . or from the pole of the stress-energy correlation function (jl.3|) (additionally 
using the equation of state to relate (|1.5jl ) must give the same result. It is much less 

1 Computation of the thermal correlation functions in the dual supergravity description was ex- 
plained in [3,4]. 
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obvious that such an agreement should persist automatically also on the supergrav- 
ity side. Thus, we regard above consistency of the hydrodynamic description of the 
black 3-branes in type IIB supergravity as a highly nontrivial check of the Maldacena 
correspondence [1] applied to near-equilibrium thermal gauge theories. 

The situation with the sound wave propagation in the hydrodynamic limit is similar 
[8] (though perhaps less dramatic compare with shear viscosity given the conformal 
invariance of the M = 4 SYM). The speed of sound can be computed from the equation 
of state as 

, dP , N 

c^ = ^, (1.6) 

where P and e are correspondingly the pressure and the energy density of the strongly 
coupled gauge theory plasma which can be extracted from the thermodynamic proper- 
ties of the black 3-branes [20]. Alternatively, all the transport coefficients can be 
read off from the dispersion relation for the pole in the sound wave channel two-point 
retarded correlation function of the stress energy tensor, for example, 

Gtt, tt (co,q) = -i [ dtdxe^ z e(t)([T tt (x), T tt (0)}) , (1.7) 



as 

Again, all these computations point to a consistent picture of a hydrodynamic descrip- 
tion of the supergravity black 3-branes 2 . 

In this paper we prove that consistent hydrodynamic description of black 3-branes 
persists even once one include leading a' correction to type IIB supergravity from string 
theory [22-25] , which translates into finite 't Hooft coupling correction on the M = 4 
SYM side of the Maldacena duality. To appreciate the nontrivial fact of the agreement 
we point to some features of a'-corrected description of the black 3-branes: 

■ including leading order a' correction, the entropy density of the black 3-branes dif- 
fers from the Bekenstein-Hawking formula which relates the latter to the area of the 
horizon [26]; 

■ the Hawking temperature of the black 3-branes as well as their equilibrium thermo- 
dynamic quantities, i.e., the entropy, energy and the free energy, receives a' correc- 
tions [26,27]; 

2 Consistency of hydrodynamic description of more complicated examples of gauge theory- 
supergravity correspondence follows from [13, 17, 18,21, 19]. 
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■ unlike the supergravity approximation [20], the radius of the S 5 of the a' corrected 
black 3-brane geometry is not constant [27]. 

We find that only properly accounting for all of the above facts one finds a consistent 
picture of the a' corrected black 3-brane hydrodynamics. Lastly, we strongly suspect 
that consistency of the hydrodynamics is sensitive to the structure of the a' corrections 
in type IIB string theory. Thus our computations can be helpful is determining exact 
structure of such corrections 3 . 

The paper is organized as follows. In the next section we discuss our computational 
approach and present the results. In section 3 we apply the general computational 
scheme to the evaluation of the dispertion relation of the shear quasinormal mode in 
black 3-brane geometry without a' corrections. This was previously discussed in [16], 
though our approach highlights the use of the effective action rather than equations 
of motion 4 . In section 4 we discuss the main computational steps leading to the shear 
and the sound wave lowest quasinormal modes dispertion relations. 



2 General computational approach and the results 

In the context of gauge theory-string theory correspondence [1] poles of the finite tem- 
perature two-point retarded correlation functions of the stress-energy tensor can be 
identified with the quasinormal frequencies of the gravitational perturbations in the 
background string theory geometry [16]. Strictly speaking, such an identification has 
been made in the supergravity approximation to gauge theory-string theory correspon- 
dence, but as it is derived from the standard prescription for the computation of the 
correlation functions [28,29,3], we expect it to be valid beyond the supergravity ap- 
proximation. In this paper we extract a'-corrected transport coefficients (jl.lj) from the 
a'-corrected dispertion relation for the lowest shear quasinormal mode (jl.4|) and the 
lowest sound wave quasinormal mode (|1.8|) in the a'-corrected black 3-brane geome- 
try [26,27]. 

We start with the tree level type IIB low-energy effective action in ten dimensions 
taking into account the leading order string corrections [22-25] 



1 



16itG 



10 



d w xy/=g 



R - \{dct>) 2 - -^(F 5 ) 2 + + 7 e~l*W + .. 



(2.1; 



3 We hope to report on this elsewhere. 

4 Using equations of motion is technically much more complicated in the presence of a' corrections. 



4 



7=^C(3)(a') 3 , 

where 

Tir /-thmnk r~t /~i rsp^iq . 1 /~ihkmnr~i /~i rsp s-^q /r> q\ 

^ _ U O pmn gO h O rsfc -+- -O O pgmn O h O rsfc . {Z.Z) 

As in [26, 14] we assume that in a chosen scheme self-dual F5 form does not receive 
order (a') 3 corrections. In ()2.1|) ellipses stand for other fields not essential for the 
present analysis. 

We represent ten dimensional background geometry describing 7-corrected black 
3-branes by the following ansatz 

, 2 



ds 2 10 =g$ dx»dx u + c\ (dS 5 )' 



2 (2-3) 
\dt 2 + 4 (dx 2 + dy 2 + dz 2 ) + c 2 3 dr 2 + c\ (dS 5 ) , 



\2 • 



where Cj = Cj(r) and (dS 5 ) is a metric on a round five-sphere of unit radius. For the 
dilaton we assume <p = <f>(r) and for the five-form 

F 5 = JF 5 + *.F 5 , T b = -4 dvoZ s « . (2.4) 

In ()2.4|) the 5-form flux is chosen in such a way that 7 = solution corresponds to 
C4 = 1. To leading order in 7, solution can be written explicitly [26,27] 

4\ 1/2 _ 



(4\ J-/ 2 5 
e-3^(l + a + 46) 

5 

c 2 =re 3 ; 



(2.5) 



r[l-# 



c 4 =e , 



where to order 0(j 2 



15rp /„_rn r, 



7 



254 - 794 + 25 



'=7^(5^-1^ + 5), (2.6) 



15rg / r 4 



v =7 1 + — 

The dilaton <fi also receives 7 corrections, <p oc 7 [26]. It is easy to see that to order 
0(7) gravitational perturbations do not mix with the dilaton perturbation; moreover 
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to study gravitational perturbations we can consistency set = 0. The Hawking 
temperature corresponding to the metric ()2.3|) is [26] 

T = T (1 + 15 7 ) = - (1 + 15 7 ) • (2.7) 

7T 

Next, consider perturbation of the five dimensional metric g$ ()2.3|) 

9$ - 9$ + V , (2-8) 
where it will be sufficient to assume that 

V = M*> ^ r ) = e ~ Mqz MO • ( 2 - 9 ) 

With the metric perturbation ansatz ()2.9|) we have 0(2) rotational symmetry in the xy 
plane. The latter symmetry guarantees that at the linearized level the following sets 
of fluctuations decouple from each other [6, 16] 

{hxy}, {hxx-hyy}, (2.10) 



h xz , h xr "\ , {J^tyi h"yz i kyr} i (2-H) 



\httj ht Z} hf T} h aa — h xx + hyy, h zz , h zr , /i rr } . (2.12) 

Scalar channel fluctuations (|2.10j) were studied in [14] leading (using the Kubo relation 
(jl.2j) ) to the following prediction for the shear viscosity to the entropy density ratio 

2=1 (1 + 1357). (2.13) 
s An 

In this paper we study shear channel ()2.11j) . and the sound channel fl2.12j) fluctuations. 
Effective action for the fluctuations (j2.11j) and ()2.12j) can be obtained by expanding 
the supergravity action (j2.1|) around the background 5 ()2.3|) to quadratic order in h^ u . 
Though we can always choose the gauge 



Hr 



5 There is a subtlety in evaluating the action with a self-dual 5-form background. The correct way 
to do this is to assume that Fs has components only along S and double that contribution in the 
lOd effective action [30]. 
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doing so on the level of the effective action for the fluctuations would lead to missing 
important constraints, i.e., equations of motion coming from the variation of the action 
with respect to {h tr , h xr , h yr , h zr , h rr }. As we explicitly demonstrate on a simple 
example in the next section these constraint equations are crucial in decoupling gauge 
invariant fluctuations. Rather, the correct way is to impose the gauge fixing (j2.14j) on 
the level of equations of motion for the fluctuations. 

Without loss of generality, for the shear channel we consider metric perturbations 
{htx, h xz , h xr }- Imposing the gauge condition h xr =0 on the equations of motion and 
further introducing 

h tx (r) = r 2 H tx (r) , h xz (r) = r 2 H xz (r) , (2.15) 
we find that the shear channel gauge invariant combination [16] 

Z shear = QH tx + UjH xz (2-16) 

decouples to order 0(7). The spectrum of quasinormal modes is determined [16] by 
imposing the incoming wave boundary condition at the horizon r — > r +, and the 
Dirichlet condition at the boundary r — > +00 in the 7-deformed black 3-brane geometry 
(J2.3j) on Z s h ear . The main steps of the computation are discussed in section 4.1. For 
the lowest shear quasinormal mode (in the hydrodynamic approximation) we find to 
order 0(7) 

to = -1 T v q 2 + 0(q 3 ) , (2.17) 



where 



r„ = \ + 6O7 + o( 7 2 ) , (2.18) 



and we additionally introduced 6 



ra = — , q (2.19) 
From IjDij) . J2ZZj), (l2~m (12131) we find 

!L = TV = Tx^ x r r , = i-(l + 135 7 + 0(7 2 )) , (2.20) 

in precise agreement with ()2.13j) reported in [14]. 

6 As will be clear from the discussion in section 4, in and q are the natural dimensionless parameters 
describing quasinormal modes. 
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There is additional subtlety in computing the lowest quasinormal frequency in the 
sound channel. Similar to the shear channel 7 we first derive from the effective action 
for the fluctuations equations of motion, and after that impose the gauge condition 



Hr 



h zr = hrr = . (2.21) 



As explained in [17], for a general five- dimensional Einstein frame background geometry 
with the metric 

dsl = -c\ dt 2 + cj (dx 2 + dy 2 + dz 2 ) + cl dr 2 , (2.22) 

with di = Cj(r), in the gauge (|2.21jl . and reparameterizing metric perturbations (|2.9J) 

as 

htt — c\ H tt , h tz = c 2 H tz , h aa = c 2 , H aa , h zz = c 2 , H zz , (2.23) 
the gauge invariant gravitational perturbation is given by 



Zsound = 4 — H t z + 2 H zz — H aa 1 + 2 n ~7n H tt , (2.24) 

U) \ uj z c 2 c 2 / uj z c\ 

and thus (in the absence of matter sector) must have decoupled equation of motion. 
In the absence of 7-corrections, the ten-dimensional Einstein frame reduces to the five- 
dimensional Einstein frame, so in defining Z soun <i we can simply take q = q. This is no 
longer the case with 7 7^ 0, as in this case the S 5 warp factor oc c 2 is no longer constant. 
Indeed, we find that defining Z sound as in (|2.23|) produces decoupled equation of motion 
only after q are rescaled as appropriate for the five- dimensional Einstein frame, namely 

c, = cT ^ . (2.25) 

Again, the spectrum of quasinormal modes is determined by imposing the incoming 
wave boundary condition at the horizon, and the Dirichlet condition at the boundary in 
the 7-deformed background geometry (|2.3j) on Z soun d- For the lowest shear quasinormal 
mode (in the hydrodynamic approximation) we find to order 0(7) 

to = c s q - i T sound q 2 + C(q 3 ) , (2.26) 

where 



r,=^ + 4o 7 + o( 7 2 ) 



(2.27) 



7 The main computational steps are presented in section 4.2. 
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Given IJDty . pTjl . (j2~T3j) . (jQU)) we conclude from (|2~2Tf) that the bulk viscosity of the 
strongly coupled M = 4 plasma is 



Of course, as finite 7-corrections translate into 't Hooft coupling corrections on the 
gauge theory side of the Maldacena correspondence, from the field theory perspective 



independent of the 't Hooft coupling. We showed here that the dual string theory 
description reproduces this fact as well, albeit in a highly nontrivial fashion which is 
moreover consistent with shear viscosity computation ()2.2()j) and alternative analysis 
in [14]. 

3 Diffusion constant of the black 3-branes hydrodynamics: 
the effective action approach 

Consider the shear channel gravitational perturbations {ht x , h xz , h xr } in the absence 
of a' corrections, i.e., setting 7 = 0. This was previously discussed in [6], where 
equations of motion for the fluctuations {h tx , h xz ] in the gauge h xr = were derived as 
perturbation of the full type IIB supergravity equations of motion. Such an approach 
becomes technically very difficult in the presence of 7 corrections: one needs to derive 
equations of motion for the deformed effective type IIB supergravity action (|2.1|) . In 
the latter case we find it much easier to derive first the effective action describing the 
fluctuations, and then derive the equations of motion from this action. The effective 
action for the fluctuations can be obtained by simply evaluating (|2.1j) to quadratic 
order in metric perturbations ()2.8|) . The important point we want to stress here is that 
the gauge fixing condition h xr = can not be imposed on the level of action. If we 
do this, we obtain two second order ODE's (coming from variation of the action with 
respect to {h tx , h xz } ) 



C = 0(7 2 ) • 



(2.28) 



(given that the latter is conformal) we immediately conclude that c 2 s 



I and C = 




) 



(3.1) 



where H... = H...(x), and all the derivatives are with respect to 




(3.2) 
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It is easy to see that given (|3.2|) equation of motion for 

Z shear {x) = qH tx (x) + XvH xz (x) (3.3) 

does not decouple. On the other hand, if we impose the gauge fixing condition h xr = 
on the level of equations of motion, we obtain an extra constraint equation coming 
from the variation of the effective action for the fluctuations with respect to h xr 

= tt>H' tx + qx 2 H' xz . (3.4) 

Notice that ()3.4|) is consistent with ()3.1|) . Given (J3.4|) we can now obtain the decoupled 
equation of motion for Z shear 

_ x 2 q 2 + xo 2 xo 2 -x 2 q 2 

~ ^shear + ^2 _ ^ 'shear + ^ _ ^3/2 ^shear ■ 

The incoming boundary condition at the horizon (x — > 0+) implies that 

^ shear i^) ^ ZghearyE) i (3-6) 

where z s h ea r(x) is regular at the horizon. Without loss of generality we can assume 

Zshear 1 ; 

the spectrum of quasinormal frequencies is then determined by imposing a Dirichlet 
condition at the boundary [16] 



Zshear 



0. (3i 



In the hydrodynamic approximation (tr> <C 1 and q < 1) the solution can be written 
in the ansatz 

Zshear = + < Zshear + 0{tf) , (3-9) 

where zfj ear , z^ ear are invariant under the scaling to — > Arc, q — > Aq with constant A. 
Substituting (|3.9j) into ()3.5|1 . and we find [6,16] 

^ (0) - 1 z {1) - -1 (3 10) 

which from (J3.8j) determines the lowest shear quasinormal frequency as [16] 

Xo = -i l -q 2 . (3.11) 
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4 Transport properties of black 3-branes at 0{{a') 2> ) order 



The general computational scheme for deriving equations of motion for the metric 
perturbation and decoupling the gauge invariant combinations for these perturbations 
is explained in section 2. The analysis is straightforward, though quite tedious. As in 
the previous section, computations are simplified using the radial coordinate x, defined 
by (jH.2j) . Both for the shear mode (|2.1fij) and the sound mode 1)2.24)1 gauge invariant 
combination of metric perturbations we find that the corresponding equations of motion 
decouple. These equations can be expanded perturbatively in 7, provided we introduce 



Z shear ^shearfl T ^shearA ) ; , . 

2 ( ^ 

Z sound Z sound ,0 T ^ sound,! 0\! ) • 

The incoming wave boundary conditions are set up at the level of the leading order in 
7, thus it is not a surprise that a natural dimensionless frequency to and a momentum q 
are introduced (see eq. 1)2.19)0 with respect to T , rather than the a'-corrected Hawking 
temperature T of the black branes (J2.7)) . 

4.1 Shear quasinormal mode 

For the shear channel fluctuations we find 



(4.2) 



x 2 q 2 + w 2 xo 2 -x 2 q 2 

-^shearfi + ^ _ ^2) 'shear ,0 + ^ _ ^2)3/2 ^hear,0 , 

_ x 2 q 2 + W 2 tv 2 -x 2 q 2 

u ^ shear,! ' _ x 2^2\ ^ shear,! ' _ ^2)3/2 ^shear,! T ^shearfl , 

where the source J s hear,o is a functional of the zero's order shear mode Z shear $ 

j _/7(4) d 4 Z s hear,0 „(3) d 3 Z s hear,0 M2) d 2 Z s hear,0 Ml) dZ s hear,0 

Jshearfl — ^ shear ^4 ^ ^ shear ^3 ^ ^ shear ^2 ^ ^ shear 

, M°) 7 

' ^shear ^shearfi ■ 

(4.3) 

The coefficients are given explicitly in appendix A. In the hydrodynamic approx- 

imation we look for the solution for Z s h ea ri m the following ansatz 



Zghearfi ~X ^ (zf h \ arfi + ^4w,0 + 
Zshear,! (^L,! + ^SL,1 + ®W 



(4.4) 
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where z^ ear i are regular at the horizon, and satisfy the following boundary conditions 



z 



(o) 

shear,0 



y (l) 

" 'shear ,0 



Z 



(0) 

shear, 1 



Z 



(1) 

shear,! 



Explicit solution of (|4.2|) subject to boundary conditions ()4.5|) takes form 



z 



(o) 

shearfi 



1- 



1 q 



Z shear,0 2 tt> ^ 



(4.5) 



(4.6) 



2; 



(o) 

shear, 1 
(1) 

shear, 1 



25 
16' 



(x 4 -4x 2 + 5) , 



32qm 



q 2 (-240 - 1565x 2 - 860x 4 + 695x 6 ) 



(4.7) 



+ 16ro 2 (594 - 264x 2 + 43x 4 )^ . 



Imposing the Dirichlet condition on x m Z s f iea r,o a t the boundary determines the lowest 
shear quasinormal frequency (|2.17|) . 



4.2 Sound wave quasinormal mode 

For the sound channel fluctuations we find 



For the sound channel fluctuations we 
n =7 „ _ (3x 2 -2)q 2 + 3« 3 2 

soundfl 2 .(_ 3ro 2 + ( x 2 + 2 ) q 2) ^undfl 

q 4 x 2 (x 2 + 2) - 2q 2 ro 2 (2x 2 + 2) - 4x 2 (l - x 2 ) 3 / 2 q 2 + 3m 4 

x 2 (1 - x 2 ) 3 / 2 ((x 2 + 2)q 2 - 3m 2 ^ s ° und '° ' 

-A„2 , q^2 



n =7 » _ (3x 2 - 2)q 2 + 3m 2 

sound,! s (_3 m 2 + ( x 2 + 2 ) q 2) '^.l 

q 4 x 2 (a; 2 + 2) - 2q 2 ro 2 (2x 2 + 2) - 4x 2 (l - x 2 ) 3 / 2 q 2 + 3m 4 

X 2 (l _ ^2)3/2(^2 + 2 ) q 2 _ 3m 2) Z — «*1 + J -«^.0 ' 

(4- 



where the source J S ound,o is a functional of the zero's order sound mode Z SO und,o 



j _r>( 4 ) d 4 Z soun dfi ^(3) d 3 Z sounc ifl -(2) d 2 Z sounc ifl dZ soun d,o 

Jsoundfi — L- sound ^4 r ^ SOU nd r ^ SOU nd ^2 ^ sound 



, /»(0) 7 



(4.9) 
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The coefficients C^ und are given explicitly in appendix B. In the hydro dynamic approx- 
imation we look for the solution for Z SOU nd % in the following ansatz 



Z S ound,0 -X im (zfo Undfl + i(\Z^md,0 + , 

Z SO und,i =x lt0 (zfc Undl + i(\zf \ ndl + 0(q 2 , 



(4.10) 



where z^ und ^ are regular at the horizon, and satisfy the following boundary conditions 



z {0) 

soundfl 



1 z {1) 

x ' ^soundfl 



= z {0) 

sound,l 

x^0 + 



= z {l) 

sound,l 

x^0 + 



x^0+ 

Explicit solution of (|4.8j) subject to boundary conditions (|4.11|) takes form 



= 0. (4.11) 

x^0 + 



(o) _ 3ro 2 + {x 2 -2)q 2 (1) 2toqx 2 



'soundfl _ 2q 2 ' soundfl 3^,2 _ 2q 2 ' 



(4.12) 



z {0) 

sound,! 



5x 2 



q 4 (2404 + 446x 2 - 4164x 4 + 2006a; 6 ) 



16(3tt> 2 -2q 2 ) 2 

- 3m 2 q 2 (1588 + 183x 2 - 2072x 4 + 1003x 6 ) + 45ro 4 (5 - Ax 2 + x 4 )^j , 

Z ^ = 8q(3J-2q 2 ) 2 ( " 13344 + 5846x2 " 4520x4 + 

- 3ro 2 q 2 (-9744 + 5035x 2 - 2604z 4 + 867x 6 ) 

- 36ro 4 (594 - 264a; 2 + 43x 4 ) j . 



(4.13) 

Imposing the Dirichlet condition on x %x °Z sound $ at the boundary determines the lowest 
sound quasinormal frequency (|2.26j) . 

Appendix 

A Coefficients of J s hear,o 

CZ, =45(1 - X*)' (A.l) 

snear -q 2 x 2 + tD 2 
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C$L = - 8 x^-q^ + ^) ( 5(1 ~ ^^i 280 ^ + 1018x2 + 216 ) " * 2 ( 72 

x (A.3) 

+ 3811a: 4 + 154a: 2 )) + 16(1 - x 2 ) 5/2 (qV - m 2 )(45ro 2 + 83qV) 



C lw = - _„ „ 2 „ 2 , ^x 2 ~5qV(l - x 2 )(1105x 6 - 2411a: 4 - 802* 2 - 216) 



8x 3 (-q 2 x 2 + to 2 ) 2 

- 4tt) 2 qV(-8792a: 6 + 1467a: 8 + 8045a; 4 + 60a; 2 - 1080) - m 4 (-13479x 4 + 3353a; 8 
+ 2930a; 2 + 5876a; 6 + 2520) + 16(1 - x 2 ) 3/2 (83a; 6 (4a; 2 + l)q 6 - m 2 qV(497a: 2 - 372) 

- q 2 roV(162a: 2 + 353) + 45tt) 6 (3a: 2 + 2))^ 

(A.4) 

CSL = - g^^ + tt 2 ) (l52q V(-l + x 2 ) - 1480m 2 qV(-l + a; 2 ) 
+ 4x 2 q 2 (-242tt) 4 + 242ra 4 a; 2 - 75a; 8 + 100a; 6 ) + 10ro 2 (-36tt) 4 + 36m 4 a; 2 + 45a: 8 

- 80a; 6 + 25a: 4 ) + (-5q 4 a; 6 (-831x 2 + 274 + 677a: 4 ) + 2ra 2 qV(-1858a; 2 + 4453a: 6 

- 5295a; 4 + 3600) - (4320 + 1250a; 2 + 6097a: 6 - 10467a: 4 ) tt) 4 )(l - x 2 )' 1/2 ^j 

(A.5) 



B Coefficients of J S0U nd, 



o 



(4) _ 2q 2 (183x 2 + 29) + 333lt) 2 

^sound 3( 2q 2( a .2 3m 2) ^ X ) \ D ' L ) 



,(3) 2(x 2 - l) 3 

Sound 3a .( q2 ( a .2 + 2 ) _ 3ra 2 )(2q 2 (a: 2 - 1) + 3ro 2 ) 2 



4q 6 (x 2 - l)(872a: 6 + 2486a; 4 



- 567a: 2 - 58) - 6ra 2 q 4 (557a; 6 - 7955a: 4 + 2386a; 2 - 106) - 9ra 4 q 2 (4237x 4 - 3875a; 2 
+ 386) - 2997tt) 6 (9a; 2 - 1)^ 

(B.2) 
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Cg2 ° und = (24x 2 (q 2 (a: 2 + 2) - 3ra 2 ))(2q 2 (a: 2 - 1) + 3m 2 ) 3 ( 8(17357xl ° + 64126x8 

- 125343a; 6 + 33528a: 4 + 4468a: 2 + 464) (x 2 - l) 3 q 8 - 12ra 2 (12419x 10 - 416676a: 8 
+ 639279a; 6 - 205814a: 4 - 18792a: 2 + 384) (x 2 - l) 2 q 6 - 18m 4 (a; 2 - 1) (183535a; 10 

- 1069144a: 8 + 1301579a: 6 - 428618a; 4 + 14112a; 2 + 3936)q 4 - 27tt 6 (a: 2 - l)(359001a: 

- 732179a: 6 + 440634a: 4 - 71592a: 2 - 4864)q 2 - 243tt 8 (a; 2 - 1) (22327a: 6 - 28997a: 4 
+ 6974a: 2 + 296) - 16(2q 2 (x 2 - 1) + 3ro 2 )(q 2 (a: 2 + 2) - 3ro 2 )(l - a: 2 ) 5/2 (4a: 2 (a: 2 - 1) 

x (125a; 2 - 87)q 6 + 4ra 2 (273x 4 - 86a; 2 + 29)q 4 - 3tt? 4 (159x 2 - 164)q 2 - 999tt? 6 )^ 



(B.3) 



C — = 24,3 (2q2( , 2 _ 1) + 3tt , 1 2 )3 (q2(x2 + 2) _3 t , 2 )3 (8(461207a; 44 + 1124342a; 42 

- 1565613a: 10 - 408216a; 8 + 1082544a; 6 - 675120a: 4 + 31312a; 2 - 1856) (a: 2 - l) 3 q 12 

- 12ro 2 (234845a; 14 - 6242508a; 12 + 5011337a; 10 + 4244306a: 8 - 6922564a: 6 
+ 3818200a; 4 - 275392a: 2 + 2176)(a; 2 - l) 2 q 10 - 18tt) 4 (x 2 - l)(3115141a; 14 

- 10599760a: 12 + 3749965a; 10 + 13178414a: 8 - 17763088a: 6 + 9000248a; 4 - 777360a: 2 

- 16960)q 8 - 27tt) 6 (-52480 + 19140582a: 8 + 10819632a: 4 - 7115740a: 12 + 3171763a: 

- 23084512a; 6 - 1929885a; 10 - 981760x 2 )q 6 - 81ra 8 (313560a: 2 + 591463a; 12 

- 1966197a: 8 - 476182a; 10 - 2942744a: 4 + 4503220a; 6 + 30880)q 4 - 243tt) 10 (-8416 
+ 96733a; 10 - 43280a; 2 + 303614a: 4 - 158080a: 8 - 215771a; 6 )q 2 - 2187m 12 (7105a: 8 

- 7980a; 6 + 296 + 426a; 2 + 1353a; 4 ) - 16(q 2 (a; 2 + 2) - 3ro 2 )(2q 2 (a: 2 - 1) + 3ro 2 ) 
x (1 - x 2 ) 3/2 (4x 2 (x 2 - l)(2863x 8 + 5423a: 6 - 8425a; 4 - 1988a: 2 - 348)q 10 

- 2ro 2 (4979a; 10 - 83715a; 8 + 87128a; 6 + 16642a: 4 + 3624a: 2 + 232)q 8 - 3m 4 (192 

+ 38586a; 8 - 110297a: 6 - 4976a: 2 - 22460x 4 )q 6 - 9tt) 6 (20610x 6 + 4957a: 4 - 148a; 2 

- 984)q 4 - 27tt) 8 (206a; 4 + 1561a; 2 + 608)q 2 + 8991ro 10 (l + 4a: 2 )) 



(B.4) 
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Cfound ~ 2 4x 4 (2q 2 (;r 2 3tt ,2)3( q 2( a .2 + 2 ) - 3 tt) 2 )3 ( G^ 643 ^ ^ 

x (a; 2 + 2) 3 q 16 + 32x 2 tt) 2 (5747a; 6 - 31591a; 4 + 25688a: 2 + 696)(a; 2 - l) 2 (a; 2 + 2) 2 q 14 

- 16(x 2 - 1) (343159a; 18 + 359756a; 16 - 2430770a: 14 - 19905tt> 4 a; 12 + 1208858a: 12 

- 216211ra 4 x 10 + 2130051a: 10 - 1974730a; 8 + 212468tt) 4 x 8 + 789188raV + 317076a; 6 

- 704624roV + 46600a: 4 - 43376roV + 928tt) 4 )q 12 + 24tt) 2 (x 2 - l)(385565a; 16 

- 4651541a: 14 + 3705393a: 12 + 6897979a: 10 - 73557ro 4 a: 10 - 7874224a; 8 + 121796ro 4 a: s 
+ 1319004a; 6 + 803308tt) 4 a; 6 - 855456tt) 4 a; 4 + 217824a: 4 - 96080tt»V + 1088ro 4 )q 10 
+ 36tt) 4 (x 2 - 1) (2216157a: 14 - 3463628a; 12 - 8055387a; 10 - 26682tt) 4 a; 8 + 12127416a: 8 

- 2103232a: 6 - 408776raV - 547176a; 4 + 579864raV + 117600raV + 8480ro 4 )q 8 
+ 54tt) 6 (26240ra 4 + 10762468a: 8 - 803584a; 4 + 2980096a: 12 + 956891a; 14 - 462808a; 6 

- 13364213a: 10 + 122304tt) 4 a; 4 - 290744ro 4 a; 6 + 58560roV + 83640tt)V)q 6 

- 162m 8 (15440m 4 + 355529a: 12 + 1905537a: 8 - 2260254a: 10 - 303468a; 4 + 349906a; 

- 2056ro 4 a; 4 - 15120tt) 4 x 6 + 1736raV)q 4 - 486m 10 (-4208m 4 + 150699a; 10 + 45576a: 

- 93108a: 8 - 111717a: 6 + 2684roV + 1524tt)V)q 2 + 4374ro 12 (-148tt) 4 - 400a; 6 
+ 148roV + 125a: 4 + 225a: 8 ) - (q 2 (x 2 + 2) - 3ro 2 )(l - a: 2 )- 1/2 (8x 4 (117035x 10 

+ 280248a: 8 - 600005a; 6 - 185858a: 4 + 519204a: 2 - 114424)(a; 2 - l) 2 q 12 - 4ro 2 (a; 2 - 1) 
x (403063a; 14 - 5091874a; 12 + 4346503a: 10 + 5410988a: 8 - 6467736a: 6 + 1273520a; 4 
+ 35760a; 2 - 7424)q 10 - 6ro 4 (10631928x 8 + 2466759a: 14 + 1592283a; 10 - 7034846a; 12 

- 8874092a: 6 + 974760a: 4 + 195888a; 2 - 1280)q 8 - 27tt) 6 (-797920a; 8 + 485056a; 4 
+ 475731a: 12 + 927244a: 6 - 918927a: 10 - 23040 - 115744x 2 )q 6 - 27ro 8 (-269131a: 8 
+ 473495a: 10 - 1348292a; 4 + 908224a; 6 + 70400 + 149104x 2 )q 4 - 81tt 10 (314448a; 4 
+ 56787a; 8 - 321287a: 6 - 26560 - 30588x 2 )q 2 - 729ro 12 (1184 + 2665a; 6 - 4027a; 4 

+ 778a; 2 )) 



.6 

A 



(B.5) 
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